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Abstract 

, We propose a simple parameterization of the two-point correlator of hadronic electromagnetic cur- 

' rents for the evaluation of the hadronic contributions to the muon anomalous magnetic moment. 

The parameterization is explicitly done in the Euclidean domain. The model function contains a 
\ phenomenological parameter which provides an infrared cutoff to guarantee the smooth behavior 

^ • of the correlator at the origin in accordance with experimental data in e~^e~ annihilation. After 



fixing a numerical value for this parameter from the leading order hadronic contribution to the 
muon anomalous magnetic moment the next-to-leading order results related to the vacuum polar- 
ization function are accurately reproduced. The properties of the four-point correlator of hadronic 
electromagnetic currents as for instance the so-called light-by-light scattering amplitude relevant 
for the calculation of the muon anomalous magnetic moment are briefly discussed. 



1 Introduction 



The description of strong interactions based on QCD proves to be very successful for processes at large 
energies where the coupling constant is small due to the property of asymptotic freedom Q. This 
makes perturbation theory (PT) computations reliable. At low energies the problem of strong coupling 
prevents using QCD as an unambiguous theoretical tool for computations of physical observables and 
various phenomenological models are introduced. These models are inspired by QCD but it is difficult 
to establish a quantitative relation between the underlying theory and a model used in practice. For 
instance, the chiral perturbation theory (ChPT) for Goldstone modes is a very general model in the 
sense that this effective theory provides the expansion of the full QCD amplitudes at low energies 
according to symmetry principles without using special assumptions. In this respect the amplitudes 
of ChPT represent the low energy limits of the exact QCD amplitudes based on current algebra and 
collect the current algebra results in a compact way [0. Of course, for the applicability of ChPT 
one should guarantee that the expansion in the energy is well convergent (at least explicitly). ChPT 
is very practical in describing interactions of pions (as lightest hadrons) with nucleons or resonances 
at low energy in the small momentum (and mass) expansion ^. Thus the description of strong 
interactions at low energies relies on phenomenological models with explicit introduction of elementary 
hadron fields or on the closely related approaches based on general principles of analyticity, unitarity 
and symmetry ||5| . A general idea of linking this approach for the description of hadrons at low energies 
with QCD is the concept of duality which means that the description of inclusive observables which 
are sensitive to the contribution of many particles is simpler than that of exclusive processes and can 
be represented by almost free fermions or weakly coupled quarks This concept works well for 
infrared (IR) soft observables in r-decays and other sum rules where the limit of massless quarks is 
nonsingular 0, ^ ^ H) 12]- For the IR sensitive observables the realization of the duality concept 
for the light modes is not quite straightforward since the IR cutoff explicitly enters the calculation. 
In such cases the IR cutoff is usually taken from experiment such as the mass of a real hadron. 

Strong interactions at low energies play an important role in precision tests of the Standard Model 
as a whole. They are involved in the evaluation of the CP-violating structure in the electroweak 
sector through hadronic matrix elements, the computation of which constitutes a main obstacle for 
the progress of a quantitative analysis of the nonleptonic kaon decays relevant for determination of 
the quark mixing matrix ||T3| , |l4| and mixing of neutral pseudoscalar states In other tests of the 
Standard Model strong interactions enter as small corrections to very accurately measured observables. 
Examples of high precision observables of the Standard Model with important contributions from 
strong interactions (as corrections to the leading leptonic effects) are the running electromagnetic 
(EM) coupling constant at the scale of the Z-boson mass and the muon anomalous magnetic moment 
(MAMM). The numerical values of these quantities put some constraints on the Standard Model 
parameters and can also serve as triggers of new physics beyond the Standard Model. 

A numerical value of the muon anomalous magnetic moment (MAMM) is measured experimentally 



with high precision |lq, 17 1. The value presented in a recent review |T8[ reads 



a^^P = 116 592 023(151) x 10"" (1) 

with an uncertainty of 151 x 10"^^. In future experiments a goal is set to reach the accuracy of 
40 X 10"^^. In theoretical computations the leading contribution to the MAMM is given by 

a 
2^ 



a?'" = ^ (2) 
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(first calculated by Schwinger) where a is the fine structure constant = 137.036... ||T^. The 
MAMM is sensitive to the IR region of integration in perturbation theory diagrams. For the leptonic 
(QED) part of the contribution this is reflected in a strong dependence on the electron mass. Theoret- 
ically the purely leptonic part is computed in perturbative QED with finite lepton masses which leads 
to a function A{me/ra^,m^/mr) which is known analytically to three loops. As mg <C <^ rrir 
the ratios me/fn^ and mi^Jnir are small and the function A{me/m^,m^/mr) can be expanded in 
these ratios to simplify calculations. The contribution of the muon leads to diagrams with a single 
scale that makes it simpler to compute. The complete analytical calculation is technically very 
complicated already at the level of three loops [^]. The nontrivial diagrams in higher orders were 
computed numerically |20, 21 1. The present value of the QED contribution to the muon anomalous 



magnetic moment reads g^, ^ (as a review see |18| ) 

^QED ^ 705.7(2.9) x 10"^^ (3) 

Computation of the electroweak (EW) corrections to the MAMM has also been performed in pertur- 
bation theory [^] . The EW contribution is now known with two-loop accuracy (as a review see jl^ ) , 



aj^w ^ 152(4) X 10-^1. (4) 

The hadronic contribution to the MAMM is sensitive to the infrared region and cannot be computed 
in perturbative QCD with light quarks. The current masses of light quarks are too small to provide 
a necessary infrared cutoff and explicit models of hadronization are required for a quantitative analy- 
25, 26, 27 1 . The hadronic contribution is the main uncertainty in the theoretical computation of 
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the MAMM in the Standard Model. Assuming the validity of the Standard Model for the description 
of elementary particle interactions 

«r = (5) 

the numerical value for the hadronic contribution to the MAMM in the Standard Model is given by 
= - "m™ - 4^ = (7165.3 ± 151|exp ± 2.9|qed ± 4|ew) x 10^". (6) 

The experimental error dominates the uncertainty. 

Since the hadronic contribution is sensitive to the details of the strong coupling regime of QCD 
at low energies and cannot be unambiguously computed in a perturbation theory framework the 
theoretical prediction for the MAMM in the Standard Model depends crucially on how this contribution 



is estimated |28|. In the absence of a reliable theoretical tool for the computation in this region one 



turns to experimental data on low-energy hadron interactions for extracting a numerical value |29|. In 
general terms the hadronic contribution to the MAMM is determined by the correlation functions of 
electromagnetic (EM) currents. Since a source for the EM current is readily available for a wide range 
of energies, one tries to extract these functions or some of their characteristics relevant for a particular 
application from experiment. Without explicit use of QCD the correction in the Standard Model 
is generated through the EM interaction ej^^^A^ with j^^^ being the hadronic part of the EM current. 
At the leading order (a^ in formal power-counting) only the two-point correlation function of the EM 
currents emerges in the analysis of hadronic contributions to the MAMM. At the next-to-leading order 
{a^) a four-point correlation function appears. These correlators are not calculable perturbatively in 
the region essential for the determination of the hadronic contributions to the MAMM. The leading 
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Figure 1: The leading order hadronic contribution to the MAMM, the shaded bubble indicates the 
hadronic two-point correlator 



contribution to the MAMM comes from the two-point correlator referred to as the hadronic part of the 
photon vacuum polarization contribution while the four-point function first emerges at the order, 
most explicitly as the light-by-light scattering amplitude. To avoid using QCD in the strong coupling 
mode one can extract the necessary contribution to the MAMM by studying these two correlation 
functions experimentally without an explicit realization of the hadronic EM current j'^J"*^ in terms of 
elementary fields. Another possibility which is close in spirit is to use phenomenological models to 
saturate these correlators with contributions of real hadrons at low energies ||2^, ^ 32|, 33]. There 
is also a possibility to use a concept of duality between hadron and quark-gluon descriptions modified 
for handling IR sensitive observables p3, EHl. In the following we discuss this last option. 



2 Hadronic contribution at leading order 

At the leading order in a the hadronic contribution is described by the correlator 



in terms of a single function Il^''^'^{q'^) of one variable q^. The contribution of H^^'^^q'^) to the muon 
anomalous magnetic moment (e.g. p7|) is given by 



'■r(L0) = 4>r(^) £ ^ A-(,)I,n (8) 



with a one-loop kernel of the form 



Kis) = I dx ^f ^"" f (9) 



Here Im n'^^'^(s) = Im \J^^^{q^)\q'i=s+io\ , m is a muon mass. 

The leading order hadronic contribution to the MAMM as depicted in Fig. |l| is represented by 
an integral over the hadron spectrum and no specific information about the function Im n'^'^'^(s) 
is necessary point-wise. However, a QCD approach based on light quark duality in the massless 
approximation is not directly applicable as the integral in Eq. (^) is IR sensitive and depends strongly 
on the threshold structure of the function li^^{q'^). In most applications the threshold structure 
is extracted from experiment. To the leading order in a the function Im n'^^*^(s) can uniquely be 
identified with data from e^e~ annihilation into hadrons. Introducing the experimental R'^^^{s) ratio 

nexD/ \ cT(e+e^ ^ hadrons) 2 /mN 
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reference 




a|l^^(LO; light) 


TTleff [MeV] 


rUg [MeV] 


Ref. 1 
Ref. 1 
Ref. 1 
Ref. 1 
Ref. 1 


28| 
36| 
38| 
39| 
29| 


7011(94) X 10-^^ 
7024(154) X 10-^^ 
6950(150) X 10-11 
6924(62) X 10^11 
6988(111) X 10-11 


6940(96) X 10-11 
6953(155) X 10-11 
6879(151) X 10-11 
6853(65) X 10-11 
6917(112) X 10-11 


200.5(1.6) 
200.3(2.5) 
201.5(2.5) 
201.9(1.1) 
200.9(1.8) 


178.7(1.4) 
178.5(2.2) 
179.6(2.2) 
180.0(1.0) 
179.0(1.6) 


optimistic 
conservative 


6961(43) X 10-11 
6979(114) X 10-11 


6892(44) X 10-11 
6908(116) X 10-11 


201.3(0.7) 
201.0(1.9) 


179.4(0.6) 
179.2(1.7) 



Table 1: Comparison of a selection of different determinations of the leading order contribution to 
in the literature. For the first five lines the first column indicates the reference, the second column the 
cited value. In the third column the leading order contributions from heavy quarks are subtracted. 
The sixth and seventh line show the optimistic and conservative estimate for the mean value. The 
fourth and fifth column list the resulting values for m-eff and rUg for model 1 and model 2, resp. 



and identifying it with the theoretical quantity 127rlm ni'^'i(s) one finds 



The contribution to the MAMM based on the representation given in Eq. (^) is well studied. Several 
recent determinations based on a thorough treatment of various sets of data are given in the second 
column of Table |l[ All values given in Table |l| are consistent within errors. We can take only a small 
selection of results because other results are again within errors but with mean values slightly different. 
For the same reason, the old historical result given in Ref. |28] could be cited as well. Table |l] gives 



two different mean values, namely an optimistic and a conservative average. The optimistic average 
uses the error estimates of the different results as weights, assuming that the results are independent. 
Here we use the method presented in Sec. 4.3 of Ref. ||3^. But because all treatments in the literature 
are mainly based on the close data sets taken for the analysis, this assumption might not be true. 
Therefore, we also give a conservative average which reads 

aJi^'i(LO) = 6979(114) x lO"". (12) 

This average is in agreement with the experimental result in Eq. (|^ within error bars. The statistical 
correlation of errors coming from the experimental value in Eq. (|l|) and the leading order hadronic 



data in Eq. (12) is supposed to be small as they are determined by different sources. Other errors are 
negligible. For the target experimental error of the MAMM at the level of 40 x lO"" @ the value 



in Eq. (12) is somewhat small and NLO hadronic corrections should be taken into account. In order 
to obtain a naive order-of-magnitude estimate for this contribution, we take the value of the leading 
order suppressed by a ~ 1/137 which gives roug hly 50 X 10-11. Writing 

a|i"1th = a;i-<i(L0) + a;i-i(NL0) 

and comparing with Eq. (|6|) one has (in units of lO-H) 

a|i^'^(NLO) = 7165±151|exp±2.9|QED±4|Ew- 6979(114)|LO 

= 186±151|exp±2.9|QED±4|EW±114|had. (13) 
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Assuming the statistical independence of the uncertainties one finds after adding them in quadratures 



aJi''<^(NLO) = (186 ± 189) x 10"^^ (14) 

which would not allow one to see higher order hadronic effects clearly. Assuming that the mean value 
of aJf^P in the planned experiment will not change and the target accuracy 40 x 10^^^ will be reached, 
one finds a numerical value for the NLO hadronic contribution of 

a|i^'^(NLO) = (186 ± 121) x 10"^^ (15) 

which makes the NLO hadronic effects noticeable at the level of two standard deviations. If the mean 
value of a^^P will change in the range of the present experimental uncertainty of 151 x 10~^^, the NLO 
hadronic effects can be more or less pronounced. 

To create a framework for the analysis of hadronic contributions at NLO based on duality argu- 
ments we rewrite the LO expression for the hadronic contribution to the MAMM given in Eq. (^) in 
a different form. From general principles a two-point correlator n(g^) as a function of the complex 
variable can have a cut along the positive semiaxis s > with a positive discontinuity |4^. This 
spectral condition plays a crucial role in the analysis of the structure of the two-point correlators and 
related observables The dispersion relation reads 

1 1'°° ds 

n(g^) = — 5-Im n(s) — subtractions. (16) 

vr Jo s - q^ 

For massive pions the experimental spectrum in e^e~ annihilation starts from 4m^ and the subtraction 
at the origin is possible as Im n'^^'^(s) = for s < 4m^. The dispersion representation with subtraction 
at the origin reads 

n>-(,') = I r jf'" (17) 

TT Jiml S S - q^ 

which implies the normalization condition n'^'^'^(O) = 0. Using Eqs. (1^,17) and Eqs. (^J9|) one can 
rewrite the LO contribution to the MAMM as an integral over Euclidean values of for I\^^^{q^)^ 

aJi^d(LO) = 47r2 [-U"^"^ {-£)] W{t)dt (18) 

with 

W{t) = ,2 - (19) 

Vt^ + 4m2i [t + 2m2 + ^/WT^rr^tj 



Such a representation is well known and is often written as a parametric integral p^, 44 1. 



The representation in Eq. (|1^) is suitable for the evaluation of the hadronic contributions to the 
MAMM by using experimental data, since it can be rewritten in terms of the hadronic cross section 



for e^e annihilation. The representation in Eq. (18) is more suitable for a theoretical study as 
perturbation theory should preferably be applied in the Euclidean domain. Integration by parts in 
Eq. (|18|) gives 



i / ^ir(.)Im n'^-'^(.) = / -^^^T^ F(t)dt, Fit) = / W{Qd^ (20) 
vr J Ami ^ Jo \ dt J Jt 
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Figure 2: The LO Euclidean weight function F{t) 



0.1 



t [GeV'] 



with 



^ y o I . . o . /To ; A 1 



2 + 2m2 + ^t'^ + AmHj + + Vt^ + 4m2t) 
The behavior of the function -F(t) is shown in Fig. for small and large t. It reads 



Fm^o = l- —+oit), 

2 m 



m 



F{t)\t^oo = ^ + Oil/t^). 



(21) 



(22) 



The surface terms of the integration by parts vanish because the integrand in Eq. (|18| ) satisfies the 
conditions |n'''''^(-t)| < CVt at small t and |n^'''^(-t)| < C't'^ at large t with some given constants 
C, C. 

A key physical quantity of the analysis is the derivative of the hadron vacuum polarization function 
dl\^^{—t)/dt which is closely related to the Adler function 



-t 



Jt ' 



(23) 



The quantity dU^'"^ {-t) / dt can be computed in perturbative QCD with massless quarks for large t, 



di 



12ttH 



1 + 



asit) 



IT 



(24) 



where eg is the charge of the quark in units of the elementary (electron) electric charge and Nc is 
the number of colors. Computation at small t in perturbation theory is not possible for light quarks 
with small masses as the theory enters the regime of strong coupling. The behavior of the function 
dU^^{-t)/dt for small t can be extracted from experiment where the lower limit of the spectrum is 
determined by the finite pion masses. This leads to a finite value for the function dH^^'^ {—t) / dt at 
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t = 0. Using the patterns of small and large t behavior of the function dll^^'^ {—t) / dt for the light 
modes, we suggest an interpolation function f{t) valid for all t in the form 



Writing 

m ^ (26, 

one has 

Pit) = In , p{0) = 0. (27) 

The analytic properties of the function p{t) are given by the cut along the positive semiaxis starting 
at s = A. The discontinuity across the cut is equal to one, 

r{s) = -lmp{-s-iO) = 9{s- A). (28) 

TT 

Thus the contribution to the MAMM contains an integral 

1(A) = / f{t)F{t)dt (29) 

JO 

which is the basic quantity for the theoretical analysis. The analytical expression for /(A) is available 
but too cumbersome to be presented here. This expression is used in numerical calculations. However, 
in order to understand the integral in Eq. (29) more deeply, in particular, to find where the integral 
in Eq. (^) is saturated or what region of integration is important, an approximation can be useful 
and worth mentioning. The constant approximation for the function f{t), 

/^PP'-(t) = const = /(O) = i (30) 

gives 

P^^^iA) = /(O) F{t)dt = /(O)— = — . (31) 

This result represents the leading term of the series expansion of /(A) for small m^. The series 
expansion of /(A) for small m? up to terms of order is given by 

1 /19 1 \ /77 \ 
/(A) = -v+^- + -lnvjv^+^- + 2lnvjv' + ... (32) 

with V = m?/ A. This series converges nicely for small values of v. 

The interpolation is only necessary for the light modes since the small t behavior is nonperturbative. 
Heavy quarks can be treated in perturbation theory as their masses are rather large [|T6| , 46 1. The 
contribution of c and b quarks to the MAMM reads [p5| 



a^^'^iLO; heavy) = 71(18) x 10"" (33) 

where we used rric = (1.6 it 0.2) GeV and nib = (4.8 it 0.2) GeV. The uncertainty mainly results from 
the uncertainty in the c quark mass. Using this result, the light mode contribution becomes 

a^^'^iLO; light) = aJi^<^(LO) - aJi^'^(LO; heavy) = 6908(116) x 10"". (34) 
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Figure 3: NLO contributions to the MAMM involving the contribution of the hadronic two-point 
correlator (left) and a lepton-hadron type (so-called double bubble) diagram (right) 



The values for the different cited values are shown in the third column of Table |l|. Using this exper- 
imental result one finds a numerical value for the IR parameter A of the interpolation function f{t). 
Writing A = 4mgg one obtains 

meff = 201.0 ± 1.9 MeV. (35) 

The individual results are shown in the fourth column of Table |^. The function r{s) of Eq. ( p8[ ) is 
depicted in Fig. |5| for rriefj = 201 MeV. This completes the quantitative description of the interpolation 
function for the two-point correlator of the light modes which can be used for the computation of the 
hadronic contributions at NLO. We cite this interpolation as model 1 in the following. 



3 Hadronic contribution at next-to-leading order 

The interpolation given by the function f{t) for the two-point correlator with the numerical value 
of the phenomenological parameter from Eq. (|35|) is now used at NLO. Two of the NLO diagrams 
involving the hadronic two-point correlator are shown in Fig. ^. The NLO contribution is an integral 
of Im n^'^'^(s) with the two-loop kernel K^'^\s), 



al^(NLO) = 47r ^K^^\s)lm U^^'^is). (36) 

The analytical expression for the kernel K^'^\s) is known |47]. 

Assuming that the IR scale of the hadronic spectrum Im n'^'^'^(s) is larger than m one can use 
an expansion of K^'^\s) in rm? / s under the integration sign in Eq. (|3^ ) to generate an expansion in 
m/Mh for the integral. The IR scale of the data is set by the explicit cutoff at y/s = 2mT^. Then one 
generates an expansion of aJ^^'^(NLO) in m? /m^ using the data. The data-based analysis for the NLO 



effects of the vacuum polarization type |48] gives 



aJi^'^(vac;NLO) = -101(6) x 10"". (37) 

In the proposed model the analysis is based on the explicit expression for the hadronic two-point 
correlator given in Eq. (p5|). In the model the hadronic scale is given by \/A = 2mcfj and the 
expansion of the kernel results in the expansion of the integral in the variable m'^/A or m? /m'^^. Note 
that convergence is not fast for the hadronic scale given by m.,^ or meff. The integral in Eq. (|36| ) can 
be rewritten in the Euclidean domain for I^^^^q^) in analogy to the LO treatment. For our purposes 
it suffices to use the expansion in m? /s. We present some contributions separately for a comparison 
with the results from Ref. 
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The vertex part of the kernel has an expansion ||4^_ 

^2 



-^ver ( 



S 



m? 


\ ( 223 


2 — 

s 




377r2 


367 


48 


~ 216 




s \llb2 48 216 \m?) 144 \m? ) ) 

( 13072841 8837r^ 10079 ( j_\ 2 ( \ 

+^1^^32000 240~ " ^600" j + ^ W"))^'"]' ^ ^ 

Generally, the terms of the expansion contain powers and logarithms of the variable m? / s. For pure 
powers one can use a generating integral representation with a polynomial P{x) 



1 dxP{x) rn'j^a^ f^^A^ ^ (39) 



m^x + s s ^ \ s I Jo 



A given polynomial P{x) restores the pure power expansion of Eq. (p^). For the logarithmic part the 
generating integral representation can be chosen with a polynomial G{x) of the form 

r^^^ = Gi(mV.) + G2(mVs)lnf4V (40) 
JO sx + \m^/ 



The polynomial G{x) generates polynomials Gi{x), G2{x) through Eq. (|40|). The mixture of pure 
powers due to the polynomial Gi{x) leads to a redefinition of the polynomial P{x) in Eq. (^). Using 
Eqs. (36,^) one finds the expression for the pure power part of the expansion 



i/ imn'^-<i(s)= / ^P(x)[-nt^-^(-m2x)] (41) 

vr Jo s Jo X 

which reduces to derivatives of n'^^'^(t) at the origin and gives the analytic part of the expansion in 
m/M/j. For the logarithmic part one finds the representation 

_/ ^^^^^1 imn'^-<i(5)= / [-n^-d(-mVrr)]G(x)dx (42) 

vr Jo s Jo 



which is sensitive to the entire Euclidean domain and gives the nonanalytic part of the expansion 
containing ln{m/Mh). This procedure can be performed up to any finite order in and the whole 
calculation can be organized in a way such that only Euclidean values of momenta are necessary for 
]-[had(-g2^^ Derivatives of the function U^'''^{t) at the origin emerging from Eq. (^) depend on the 
hadronic scale Mh while the muon mass enters polynomially. The integral in Eq. (|4^ ) depends on 
both the hadronic scale and the muon mass. The behavior of Il^^'^{q'^) in the Euclidean domain is 
smooth and perturbative at large momenta. The region near the origin = is a nonperturbative 
one. Thus, the basic objects that emerge in the analysis are derivatives of the function n'^'*'^(t) at the 
origin and integrals of the form 

/JV-(-«)l„'(J,). (43) 

One can use this technique to avoid any reference to the physical region. 

In model 1 given by Eqs. (2^,27) the explicit interpolation function p{t) is given in the whole 



complex t-plane. Therefore, it makes no difference how one computes the necessary integrals either in 
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the Minkowskian or in the EucHdean domain. At the formal mathematical level the calculation with 
an explicit function can be performed in the spectral representation alone. The Euclidean approach 
will lead to the same formal results. For the calculation of the basic objects emerging in the spectral 
representation 

•^..(A) = A"/;ji,l,/(^) (44) 



the recurrence relation 

Mn,p{^) = ^ lnP(A/m^) + ^Mn,p-i{A) (45) 



n n 

can be used to decrease a power of the logarithm in the integrand of Eq. (^^. The results for the first 
two powers of the logarithm are 



A^„,o(A) = -, Mn,i{A) = - ln(A/m2) + — . (46) 



-, 7W„,i(A) = iln(A/m2) + i 
n n n 

Writing the analytical expression for the NLO vertex contribution to the MAMM in the form 

3 



a;^°'^(ver; NLO; analyt) = ( - ) V{m^/A) (47) 



a 



vr 



(the index "mod" stands for model 1, "ver" represents the vertex contribution in Fig. ^ left) one finds 

, V [2,77 2 23, \ v"^ /23647 STvr^ 677, 19 , o \ . 2^ x 

Viv) = + — \nv] + — \ + \nv + — \n^v \ + oiv^) 48 

^ ' 9V18 6 / 9\^ 1152 16 144 48 / 

where o(f^) is any function that satisfies lim^j^o o(f^)/w^ = 0. For brevity we have explicitly presented 
only two terms of the expansion of the function V{m? / 1^) at small m? / A. resulting from the correspond- 
ing expansion of the kernel in Eq. (|38|). The result of evaluating the n, d, s light mode contribution 
obtained from Eqs. (47,48), together with the appropriate QCD group factor Nc{e^ + + e^) = 2 
reads 

a;f°'^(ver; NLO; light; analyt) = -191 x 10"". (49) 
The numerical integration of the kernel given in Eq. (|38|) results in the value 



a^^^°'^(ver; NLO; light) = -190(2) x 10"", (50) 



to be compared with the corresponding results from the model of Ref. p5[| which is based on the use 
of free massive fermions for computation of the interpolation function, 

a;;i^°^(ver; NLO; light) = -188 x 10"". (51) 

The total contribution of the vertex type including the heavy quarks reads 

0;^°"^ (ver; NLO) = -194(3) x 10"", (52) 

to be compared with the one obtained in the data-based approach 



Ref. [48](^g^.NLQ) ^ _2ii(5) X 10-". (53) 



Note that the leading order contribution used in Ref. [48| is different from the value in Eq. (p^). 
Therefore the direct comparison should be made with some caution. In fact, the value used for the 
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contribution 


model 1 


model 2 


model 3 


a-°<i(LO; light) 


6908(114) X 10"^^ 


6908(114) X 10-11 


6578 X 10-11 


aJf°<i(LO; heavy) 


71(18) X 10-11 


aJf°'i(ver;NLO;hght) 


-189.5(2.4) X 10-11 


-188.2(2.4) X 10-11 


-183.2 X 10-11 


a;f°'^(ver;NLO; heavy) 


-4.3(0.9) X 10-11 


a;j^°'i(db;NLO;nght) 


105.0(1.7) X 10-11 


104.9(1.7) X 10-11 


100.0 X 10-11 


a;f°^(db;NLO; heavy) 


1.1(0.3) X 10-11 


a;f°'i(vac;NLO) 


-87.7(0.7) X 10-11 


-86.5(0.7) X 10-11 


-86.4 X 10-11 



Table 2: The different vacuum polarization type LO and NLO contributions for model 1 (second 
column), model 2 (third column), and model 3 (fourth column). The values for the effective masses 
for model 1 and model 2 are taken as rries = (201. Oil. 9) MeV and rriq = (179.4ib0.6) MeV respectively, 
the cited uncertainties are a consequence of the uncertainties of these effective masses and of the masses 
of the heavy quarks. The error estimates due to these two sources can be added quadratically. 



LO contribution in Ref. |48| is larger than the value in Eq. (|3j) which would result in a smaller value 
of rUeS and consequently the larger value of the NLO vertex correction in the model calculation given 



in Eqs. (49,^). Nevertheless, the difference is within the error bars. 

For a mixed contribution of the lepton-hadron type (so-called double bubble ("db") diagram) 
shown on the right hand side in Fig. |3| we find 

a^'"^{dh; NLO; lept + had) = 106(2) x lO-n (54) 

while the data-based estimate reads 



Ref. [48] ^LQ. igp^ ^ Y^^^^ ^ ^Q7(2) x lO-n . (55) 



Thus, the results obtained with the interpolation model given in Eqs. ( p5| , p7| ) reproduce those obtained 
in the data-based approach. 

The result for the total NLO hadronic contribution of the vacuum polarization type (two-point 
correlator only) is 

a;f°'l(vac;NLO) = -86.5(0.7) x lO-n. (56) 

which has to be compared with the result shown in Eq. (^). All contributions, including the leading 
order contributions which reproduce the input of the determination of rricfj are listed in the second 
column of Table |2|. 



The model given in Eqs. (25,^) reproduces rather accurately the results for the NLO hadronic 
contributions found in the data-based analysis for the graphs related to vacuum polarization. This 
has been expected as these results are obtained by the integration of the two-point function with the 
NLO kernel K^'^\s) which has a structure very close to that of the leading order kernel K{s). These 
two functions are shown in Fig. 0. 



4 Alternatives for the model spectral function 

The interpolation function for the two-point correlator of hadronic EM currents in Eq. ( |25| ) is very 
simple. One can use more sophisticated interpolations. A formal criterion for the choice of the inter- 



12 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

s [GeV'] 

Figure 4: The LO and NLO kernels K{s) and K^'^\s) 



polation is its consistency with general principles of quantum field theory (analyticity and unitarity in 
this case). A practical criterion is its simplicity such that analytical calculations become technically 
feasible. One can turn to free field models in a search for mathematical functions that can be used in 
the interpolation procedure. For instance, the scalar or fermionic polarization functions with masses 
as free parameters can be taken as suitable candidates. The fermionic interpolation function was 
considered in detail in Ref. [35|. It is given by the expression 

nit, nig) = (^^-1^ ip{z) - 
if{z) = artanh(V^) - 1, z= ^J^^^ - (57) 



1 

The discontinuity across the cut {^m'^, oo) at t = —s — iO is given by the fermionic spectral density of 
the form 



^ , 4m^ / 2m^\ 

A pictorial representation of Pq{s) is shown in Fig. ^. The two functions f{t, mcfr)/3 and —d7T{t, mq)/dt 
coincide within 1% accuracy in the interval t = (0, m^) if the effective parameters are related through 
fncs/iriq = \/5/2 1.12. This is expected from duality considerations based on the shape of the 



discontinuity across the cut given in Eqs. (2^,5^. Since the integral in Eq. ( p9[ ) is saturated at 
the scale of the order of the muon mass \/t ~ and < meff, the closeness of the integrated 
interpolation functions implies a rather accurate equality of the resulting integrals. The generalization 
to interpolations based on scalar field theories is straightforward. The interpolation given by the 
spectral density in Eq. ( |58D will be referred to as model 2 in the following. 

At this stage it seems that hadrons have completely disappeared from the analysis and some 
artificial functions are rather arbitrarily used to compute the relevant integrals. The link to physics 
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is that for such kind of inclusive observables as the MAMM which are sensitive to a contribution of 
many hadrons the analysis can completely be done in the Euclidean domain with the only IR sensitive 
contribution coming from the region near the origin. In the Euclidean domain the contribution of 
all hadrons to the MAMM is smeared to the extent which is determined by the distance between 
the integration region and the nearest physical singularity. The hadronic singularity is taken from 
experiment as a two-pion cut with the kinematical constraint s > Am^. The integration, in fact, 
is sensitive to scales ^/i ~ 2m^. Therefore, a correctly normalized function (as for duality at large 
energies) with the appropriate IR behavior in the interval < t < 4m^ reproduces the data with 
reasonable accuracy. The IR behavior is mainly determined by the experimental scale 4m^. 

To illustrate this statement let us consider a more realistic interpolation for the vacuum polarization 
function in the Euclidean domain from the point of view of experimental hadron physics than those 
described so far. The model hadron spectrum for light modes can be chosen in the following simple 
form 

p^^{s) = 2ml6{s - ml) + 9{s - 2ml). (^9) 

This is a one-scale no-parameter model that satisfies the duality constraints from the operator product 
expansion [0, ^ The hadronic scale of the model is given by the p- meson mass mp which is 

eventually fixed from experiment |16]. We neglect small violations of fiavor symmetry for the u,d,s 
light modes and consider them to be degenerate. The spectrum from Eq. (|5^ ) gives an interpolation 
function of the form 

V, A 1 
/had(i) = fL^ + 60) 

^ ^ {t + m2)2 t + 2mj ^ ' 

The value of the interpolation function f^^^{t) at the origin t = reads 



5_ 



/had(o) = __. (61) 

Imi 



It should be compared with the value from Eq. (pO|). For mp = 769.9 MeV one finds f^^{'d) = 
4.22 GeV-2 while the data give /(O) = 1/A = 6.25 GeV 2. This is fairly reasonable given the 
simplicity of the model but is not accurate enough. Note that there is no single free parameter in the 
model given in Eq. (|5^) and the shape of the spectrum is fixed from the duality constraint at large 
energies. The weight function F{t) from Eqs. ( pO| , ^ ) determining the MAMM integral can resolve 
the behavior of the interpolation function f^^{t) for the hadron correlator at the scales of order 
mp = 105.66 MeV. However, the approximation of infinitely narrow resonance in Eq. (|59| ) is too rough 
for computing such an integral. Thus, the low energy behavior of the spectrum is not precise enough 
which leads to an insufficient accuracy of the interpolation function at small t. Formally this is seen 
in the absence of the scale 4m^ which is known to be important for the evaluation of the integral in 
the data-based analysis. Therefore the spectrum from Eq. ( ^9| ) should be corrected for this particular 
application. A dominant role in the data-based analysis is played by two-pion states. The p meson is 
a resonance in the two-pion system therefore its contribution eff^ectively takes into account the pion 
singularity as well. However, the zero width approximation is not good enough for computing the 
particular integral in Eq. (pO|). A natural modification of the spectrum is to introduce a finite width 



for the p meson. This is achieved by replacing the function 5{s — mV) by the Breit-Wigner function 



for the resonance part of the spectrum in Eq. (|5 

>ad. X _ Tpmp 
Pr [s) - 



7T (s-m2 + r2/4)2+r2m2' 
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Figure 5: s-dependence of the spectral functions pi{s) = r{s) of model 1 in Eq. (|2^ ) and P2{s) = 3pq{s) 
of model 2 in Eq. (^8|) (upper diagram), as compared to the spectral function psis) = p^^'^{s) for 
model 3 in Eqs. ( ^|6^ ). We use mes = 201 MeV, niq = 179 MeV, and the central values nip = 



769.9 MeV and Tp = 150.2 MeV [|16|. 
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p^'^^{s) = e{s-imi)9i2mf,-s)p^^''{s) + 9{s-2mp. (62) 

The interpolation function based on this spectrum will be called model 3 in the following. Fig. |5| 
shows the s-dependence of /Op'^'^(s). The expression for the resonance part of the spectrum reduces 
to 2m?p6{s — m?p) in the limit T p 0. Using the Breit-Wigner form of the spectrum for the region 
4m^ < s < 2mp one finds the contribution of the resonance to the interpolation function in the 
Euclidean domain 

2ml 

frit) = J ^1^. (63) 

The interpolation function in the Euclidean domain for the spectrum with nonzero width reads 

frit) = /r (t) + (64) 



Computing the value of the interpolation function at the origin for Tp = 150.2 MeV [16| one finds 



/had(o) = /^^<i(0) + ^ = 5.15 + 0.84 = 6.0 GeV^^ (65) 
2m^ 

instead of the result (|6^) obtained in the infinitely narrow resonance approximation. The number 
from Eq. (|6^) differs from the data-based estimate /(O) = 1/A = 6.25 GeV~^ by 4% only. Note that 
no free parameters have been used so far for the description of the hadron spectrum. The integrals 
entering the MAMM as in Eq. (p9|) are also rather close. One finds 



to be compared with 

/had ^ 0.0155 + 0.0030 = 0.0185 (67) 

where the first term in the sum is given by the resonance and the second by the continuum contribution. 
Thus, the simple and parameter-free model from Eqs. (^J6^) already gives a reasonable precision of 
about 5% for the LO hadronic contribution to the MAMM. In this sense it successfully incorporates 
experimental information necessary for the MAMM computation. 

The spectral functions taken for model 1 from Eq. (^), for model 2 from Eq. (|58|), and for model 3 
from Eq. ( |62| ) are shown in Fig. |5|in order to allow one to compare these models. Neither model 1 nor 
model 2 has a discontinuity across the positive semiaxis of the s-plane resembling the experimental 
spectrum. However, both models result in integrals over the spectrum for the respective kernels which 
are very close to the result obtained in the experimentally inspired model 3 and, eventually, to the data 
represented in Eq. (^). From the purely mathematical point of view this is related to the fact that 
the procedure of analytic continuation is an incorrectly posed problem: small variations of functions 
in the Euclidean domain can produce big variations on the cut. The t-dependence of the Euclidean 
representation by the functions f{t), —3d7r{t)/dt, and /p'^'^(i) is shown in Fig. ^. A phenomenological 
interpretation of the situation consists in a duality between hadrons and free light fermions with QCD 
quantum numbers as for the particular application related to the MAMM computation. Of course, the 
main objective of the calculation of the hadronic contribution at the leading order from experiment 
is to reach a high precision. The use of direct data seems to be superior to a parameterization of 
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Figure 6: The functions fi{t) for the three different models where fi{t) = f{t) is given by Eq 
/2(t) = -3dTr{t)/dt is given by Eq. and hit) = /^""^(t) is given by Eq. We take A = 

where meff = 200 MeV is used for the upper diagram and meff = 205 MeV for the lower diagram. The 
parameter used for 7r(t) is connected to meg by mg = 2meff/\/5. The values nip = 769.9 MeV and 

p _ 1 r;n OA/r^^r -.lor,^ /had 



150.2 MeV used in ff^""^ are taken from Ref. [|T6|. 
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Figure 7: The light-by-light contribution (left) and the two-photon Green function (right) to the 
MAMM at NLO 



the spectrum from indirect observations. However, as soon as the integral over the data is computed 
a smooth interpolation function of a simple form can be introduced in the Euclidean domain to be 
used in higher order calculations. Because this interpolation function is explicit and complies with 
the general properties of analyticity and unitarity one can find its discontinuity across the positive 
semiaxis and perform further calculations in the spectral representation as well. The analysis of NLO 
contributions along these lines shows that the data-based results are accurately reproduced |35]. 



5 Discussion of the light-by-light contribution 

Still the class of interpolation models based on field theories contains more than just a useful set of 
functions with suitable properties. It can be used, with some caution, for extrapolations to higher order 
correlation functions as well. Indeed, the interpolation models Eqs. (p7|, 57) suffice for the calculation 
of the polarization-type hadronic contributions to the MAMM related to the two-point correlator 
of hadronic EM currents. For the whole computation at NLO one needs a four-point correlator of 
hadronic currents. It appears in two instances: as a contribution to the light-by-light graph and 
a two-photon Green function (see Fig. |^. For the contributions to the two-photon Green function 
the relevant projection of the four-point correlator depends on one external momentum and the above 
discussion is applicable theoretically as the interpolation function of one complex variable is necessary. 
There is no experimental data to fix the IR scale though. For the light-by-light graph the projection 
of the four-point correlator relevant for the MAMM calculation is a function of two independent four- 
momenta ki, k2- The scalar form factors g^ikf, /c^, (fci — ^2)^) of the projection of the tensor four-point 
correlator relevant for the MAMM computation are given by 



= ^ T-,;,(A;i, k2) gl {kj, kl {ki - k^f) . (68) 



The functions g\ are functions of the three scalar variables kf, /c2, and {ki — ^2) . Here T^J^^iki^h 



2 



are tensor structures which are polynomials in the four-momenta ki, k2 and the metric tensor g^^. The 
form factors g\ depend on the IR hadronic scale which is related to the experimental masses of light 
hadrons (pions). It may effectively differ from the scale emerging in the two-point correlator. These 
IR scales cannot be found theoretically in QCD as the regime of strong coupling is not treatable 
in perturbation theory. They can perhaps be estimated in some nonperturbative approach as, for 
example, the lattice approximation. The scalar form factors g\{ki, k2, {ki — /c2)^;M/j) are integrated 
over kl, /c2 with the weight functions w{ki , k2,m) to give a contribution to the MAMM in a full analogy 
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with the two-point correlator 

gi(^klkl{ki-k2f;Mh)w\h,k2,m)d*kid^k2. (69) 



E 



The weight functions w^(ki, k2,m) are generated by perturbative diagrams and depend on integration 
momenta ki, k2 in the loops and on the muon mass m (more precisely, they depend also on the muon 
momentum p which is taken on the muon mass-shell = m?). To the leading order of expansion 
in m/Mh the integration in explicit models reduces to vacuum three loop bubbles which can be done 
analytically Q. 

The calculation of the integral in Eq. (|69|) within an hadronization procedure requires one to 
establish the analytic properties of a given form factor g\{zi, Z2,z^] M^) as a function of three complex 
variables Zi . Also one has to saturate the singularities with the contributions of hadrons in full analogy 
with the two-point correlator. Then the integral in Eq. (|69|) is calculated through the discontinuities 
on the singularities. In other words, one can use the dispersion representation for the form factors 
g\{ki, k2, [ki — k2)'^;Mh) with some spectral density and then integrate over the momenta ki, k2 
explicitly. At this point one is left with an integral over the whole physical spectrum for the form 
factors with kernels obtained after such integration. Note that this is close to the way how the actual 
integration of the light-by-light diagram was done analytically for leptons where the physical spectrum 
can be computed in perturbation theory []l9| . 

While the physical spectrum at low energies cannot be computed in QCD perturbation theory 
point-wise, the integrals over the hadronic spectrum for the three-point functions are quite accu- 
rately calculated in the sum rule approach based on duality between hadron and quark-gluon con- 
tributions [^]. As one needs only integrals of the functions glikf, k2, {ki — k2)'^]Mh) one can avoid 
constructing these functions point-wise. Assuming general requirements of smoothness one can use 
an interpolation function in the Euclidean domain. Recently it has been argued that the neutral pion 
contribution to the amplitude Eq. (^) does not give a singular dominant contribution by itself in 
the kinematical region relevant for the MAMM calculation and can be treated within the duality ap- 
proach (see also Ref. ||52|). This implies that the integration over the Euclidean domain properly 
takes into account singularities of the (74-functions related both to the pion contribution and two- 
particle cuts given by the two-pion states. Thereby all contributions at low energies of the integrand 
(spectrum) are taken into account, for example, the p meson. The integral of the spectrum can there- 
fore be approximated by an integral of some interpolation function with an appropriate IR parameter. 
As the IR scale is set by the pion mass rriTj one would expect a numerical value of the IR parameter to 
be of the same order of magnitude. In contrast to the two-point correlator Il^'^'^{q'^) the 54-functions 
related to the four-point correlator in Eq. ( ]68| ) are complicated functions of three complex variables 
which makes it difficult to find an appropriate candidate for an interpolation function obeying all 
general requirements of field theory. A simple way to find such an appropriate interpolation is to use 
a field theory for constructing the spectrum. The use of a free fermionic theory is just the way one 
can generate the interpolation function having the necessary properties. Indeed, it does not violate 
any general principles (gauge invariance, discrete symmetries, analyticity, unitarity, . . . ) as the finite 
order perturbative field theory is about the only realistic field theory that is known to obey the general 
principles. One can consider such an interpolation as an efficient adaptive integration procedure like 
the spline technique or Monte Carlo routine VEGAS [53|. It is not a full-scale approximation of QCD 



at low energies that would allow one to compute the exclusive characteristics of hadrons, i.e. to com- 
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pute the functions gl{kl, fcl, (^i — ^2)^; M^) point-wise. It is a way to find integrals of the special type 
in Eq. (^) using duality between quarks and hadrons. Note that the duality idea between s-channel 
resonances and Regge trajectories was very useful in description of hadron scattering before the in- 
vention of QCD. Thus, many model interpolations are possible in the leading order and give rather 
accurate results. However, in the higher orders of EM perturbation theory there are strict constraints 
on the interpolation functions that one can use. Indeed, it is difficult to generalize the models given 



in Eqs. (P5|,S3) to higher orders of perturbation theory. The first model is rather artificial indeed 
and is not literally realized in any field theory while the second requires to work out the full-scale 
hadronization of QCD for the four-point correlator already at NLO. This is technically difficult as the 



data-based approach shows [p3| , p4| ]. At the same time the model given in Eq. (57) for the two-point 
correlator can immediately be extended to any order of perturbation theory in the electromagnetic 
interaction as it is a free fermion theory with ruq being an IR regulator. One also knows that this 
model is accurate for the LO hadronic contributions with the IR scale taken from the data. The IR 
parameter effectively estimates the distance from the Euclidean domain of integration to the physical 
singularities of the hadronic correlation functions. The key physical point of the quantitative analysis 
of the hadronic contributions to the MAMM within this approach is that the same parameter rries, 
or niq, enters both two-point and four-point correlators. This is, of course, an assumption which is 
based on the observation that the numerical value of this parameter is close to the pion mass. The 
absence of the neutral pion pole in the g\ functions is essential for such an assumption to be valid. 
The result of the analysis of the total next-to-leading order hadronic contributions to the MAMM 
based on the fermionic interpolation with the same value of the effective IR scale Mh = rUq for the 
two-point (polarization type) and four-point (light-by- light) correlators [35| 



a;;;'°'^(NLO) = (85 ± 20) x 10"^^ (70) 



agrees with the experimental value from Eq. (15). Note that this result includes the explicit con- 



tribution of the EM correction to the polarization function at the NLO, i.e. the basic normalization 
quantity obtained from the data as given in Eq. (|3^) is supposed not to contain this type of contribu- 
tions. This can be a rather clumsy arrangement from the experimental point of view but it is more 
definite concerning the theoretical quantities. Numerically the difference is well within error bars for 
the LO contribution though. 

If the scalar field theory of charged pions given by the Lagrangian 

= \Df,Trf - mil?, D^ = d^- ieA^ (71) 

is used to generate the four-point correlator of hadronic currents and thereby the form factors g\ 
that enter as integrands in Eq. ( |69D then the IR scale is explicitly identified with the pion mass and 
can be taken from experiment. In this case one has to add the contribution of higher resonances to 
satisfy the duality constraints at large energies which would make this approach equivalent to the one 
based on explicit hadronization. The free fermionic theory with the QCD arrangement of quantum 
numbers is exceptional in this sense as it automatically complies with the duality constraints at large 
energies. One, of course, should remember that this is a model and its large energy behavior is 
accurate only up to higher order QCD corrections. For the applications of interest this is inessential. 
A technical advantage of the fermionic theory is that the analytical results for the MAMM at NLO are 
available. For the scalar theory numerical results are available at present although the calculations can 
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in principle be done analytically as well since all necessary master integrals have been found |19, p5| ]. 
They have already been used in three loop calculations. According to the hadronization picture the 
contributions of the fermionic interpolation functions at low energies should be substituted by the 
pionic ones. In the pure fermionic model with a small effective mass the replacement of the hadronic 
contributions by the model ones is effectively done at rather low energies which makes the separate 
contribution of pions small or even vanishing. 

6 Summary and conclusions 

A parameterization of the photon vacuum polarization function related to the light hadronic modes 
is described in the Euclidean domain. The model contains a single parameter which is fixed from the 
experimental result for the LO hadronic contribution to the MAMM. The model describes the NLO 
hadronic contributions of the vacuum polarization- type in agreement with existing estimates. The 
calculation of the total NLO hadronic contribution to the MAMM in a closely related model based on 
the fermionic interpolation of correlators of hadronic electromagnetic currents is also discussed. 
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Note added 



Recently the authors of Ref. |22| have updated their result concerning the neutral pion contribution 
to the light-by-light diagram. While the magnitude of the numerical value remains the same, the sign 
of the contribution has been changed [E^] . This change brings the explicit hadron-based result for the 



light-by- light contribution rather close to that obtained within the duality approach of Ref. |35| which 
has further been investigated and developed in the present paper (originally posted at the hep-ph 
ArXiv as hep-ph /0111206| ). Thus, the present theoretical result for the muon anomalous magnetic 



moment agrees with the current experimental value. 
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